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ABSTRACT 

We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 
2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemen¬ 
tal spectroscopy and polarization data from the Steward Observatory. We find that the emission lines 
are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability 
over the 5 years. This is consistent with photoionization primarily by an accretion disk, allowing us 
to use the lines as a probe of disk activity. Comparing optical emission line flux with Fermi y-ray flux 
and optical polarized flux, we investigate whether relativistic jet variability is related to the accretion 
flow. In general, we see no such dependence, suggesting the jet variability is likely caused by internal 
processes like turbulence or shock acceleration rather than a variable accretion rate. However, three 
sources showed statistically significant emission line flares in close temporal proximity to very large 
Fermi y-ray flares. While we do not have sufficient emission line data to quantitatively assess their 
correlation with the y-ray flux, it appears that in some cases, the jet might provide additional pho- 
toionizing flux to the broad line region, which implies some y-rays are produced within the broad line 
region, at least for these large flares. 

Subject headings: BL Lacertae objects: —galaxies: active — galaxies: jets — techniques: spectroscopic 
quasars: emission lines 


1. INTRODUCTION 

Blazars are radio-loud Active Galactic Nuclei (AGN) 
whose relativistic jet are pointed at sma ll angles 
with respect to our line of sight (lAntonuccil 119931 
lUrrv fe Padovanilll995lf . This orientation makes blazars 
an ideal laboratory for the study of jet physics, due 
to Doppler beaming that greatly increases the source 
brightness and decreases the variability timescale. 

Blazars are observable across the electromagnetic spec¬ 
trum. Their broad band spectral energy distribution 
(SED) has two characteristic peaks: one at low frequen¬ 
cies (infrared - soft X-ray) due to synchrotron radiation 
and one at high frequencies (MeV - TeV), likely due to 
inverse Compton scattering. 

Since the launch of Fermi in 2008, the GeV be¬ 
havior of blazars has been studied in unprecedented 
detail. Time variability studies have constrained 
emission models by identifying correlations between 
the two spectral peaks, using high cadence obser¬ 
vations with Fermi Large Area Telescope (LAT) 
and coordinated multiwavelength cam p aigns from 
radio to TeV (e.s. IBonning et al 120091 lAbdo et al~ 
120091 iD’Ammando et all 12009 : IGhisellini et al. 


1 Chancellor’s Faculty Fellow, Syracuse University, Depart¬ 
ment of Physics, Syracuse, NY 13214;jcisler@syr.edu 

2 Department of Astronomy, Yale University, PO Box 208101, 
New Haven, CT 06520-8101 

3 Department of Physics and Yale Center for Astronomy and 
Astrophysics, Yale University, PO Box 208121, New Haven, CT 
06520-8120 

4 Yale-NUS College, 6 College Avenue, East, Singapore, 
138614 

5 Steward Observatory, University of Arizona, 933 N. Cherry 
Avenue, Tuscon, AZ 85721 


et al] 

21 ) 11 1 : 

et al] 

iia 

eWdJ 

2013 : 

et al.l 


However, the thermal components of AGN, namely the 
accretion disk and broad line region, are still energeti¬ 
cally relevant. For example, the disk could, via magnetic 
interactions, co ntribute to the initia l co llimation of the 
relativistic jet dBlandford fc Pavnei I1982H . although the 
exact mechanism is not well understood. Furthermore, 
in some blazars, the accretion disk has been shown to 
contribute a significant fracti on of the radiation e nergy 
density on sub-pc scales fe.g. IGhisellini et aLll2009f ). 

While many studies of the total flux variability of 
blazars have been undertaken, similar high-quality, 
multi-epoch spectroscopic studies have been more chal¬ 
lenging. Thus, in this work we have focused our analy¬ 
sis on the emission line variability, as the emission lines 
could be an appropriate proxy for the disk emission, 
which is often swamped by jet continuum in high flaring 
states. We have another manuscript in preparation that 
will address the continuum variability of these blazars 
and what, if any, relationship can be drawn to the emis¬ 
sion line variability properties discussed here (Isler et al. 
2015, in prep). 

Early spectroscopic studies of blazars showed emis- 
sion line variabilit y on month to year t imescales (e.g., 
Ulrich et al.1 119931 iFalomo et al.1 119941 iKoratkar et al'l 
199811 . but were not carried out in conjunction with 
7-ray observations, so the impact of the jet on these 
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sources could not be easily investigated. In principle, 
the relativistic jet could contribute additional photoion- 
izing flux to the emission lines, causing significant jet¬ 
cooling within the broad line region, provided the pho- 
toionizing emission arises on smaller spatial scales than 
the broad line gas. This geometry is required by the for¬ 
ward beaming of the jet emission, very little of which 
is directed backwards. We are now able to compare di¬ 
rectly (and simultaneously) the multi-epoch optical spec¬ 
troscopic observations of the broad line region flux to the 
jet flux using Fermi , SMARTS and Steward Observatory 
data. If a relationship is found, the spatial scale of the 
7 -emitting region can be tightly constrained. 

Simultaneous emission line variability studies in 
Fermi-monitored blazars have generated mixed results. 
Among a set of similar 7 -ray and optically bright, 
variable quasars, no emission line variabilit y was de¬ 
tected in PKS 1222+216 or 4C 38 . 51 dSmith et alJ 
1201 lb iFarina et alJ 120121 : iRaiteri et alJ l2012fh By con- 
trast, 3C 454.3 had factors of two emission line varia- 
tion roughly coincident with high y-ray e mission levels 
(llsler et al.ll2013fc iLeon-Tavares et al.!l2013 1. These stud¬ 
ies underscore the importance of monitoring blazar emis¬ 
sion line variability, but represent too limited a sample of 
the total population to draw statistical conclusions about 
blazars as a class. 

We measure the emission line behavior of 7 blazars 
to investigate the presence of short timescale (weeks to 
months) emission line variability, and assess if that vari¬ 
ability is temporally related to Fermi 7 -ray flares. This 
study could provide a direct test of the contribution of 
photoionizing flux from the jet to the broad line region, 
and if a correlation is found, observationally constrain 
the location of the 7 -emitting region to be within the 
broad line region for those flares. 

In Section O we discuss the sample selection, obser¬ 
vational program and data analysis. In Section [3l we 
present the Fermi 7 -ray and emission line flux light 
curves and we define empirical line flares as well as the 
measures for statistical variability. In Section [4] we ana¬ 
lyze the optical linear polarization at the time of obser¬ 
vation as an additional measure of the non-thcrmal jet 
contribution to the optical waveband. We discuss the 
emission line behavior of the sample and its implications 
for current jet dissipation models in Section [5j main con¬ 
clusions are summarized in Section [ 6 ] The following cos¬ 
mological parameters were used throughout this work: 
H 0 = 71 km s -1 Mpc -1 , f2 m = 0.27, Ao = 0.73, and qo 
= - 0 . 6 . 

2. SAMPLE SELECTION AND DATA ANALYSIS 

Since 2008, we have carried out optical spectroscopic 
monitoring of approximately 30 Fermi 7 -ray bright 
blazars at the queue-scheduled Small and Medium Aper¬ 
ture Research Telescope System ( SMARTS ) in Cerro 
Tololo, Chile. 

2 . 1 . Sample Selection 

The SMARTS Multi-epoch Optical Spectroscopy Atlas 
(SaMOSA) was based on the original Fermi -LAT ‘bright 
source list’ released just before launch in 2008, includ¬ 
ing those blazars with declination < 20 °, given the lo¬ 
cation of the SMARTS telescopes. The original list in¬ 
cluded Flat Spectrum Radio Quasars (FSRQs) and BL 


Lac objects (BLLs), which are distinguished by whether 
broad emission lines are present at levels greater or less 
than 5 A, respectively ijAnerel fc Stockmanlll980f l. In sub¬ 
sequent years, the SaMOSA list was expanded to in¬ 
clude newly-flaring Fermi blazars, defined as having F T 
(E > 100 MeV) > lxlO -6 ph s -1 cm -2 , and flaring FS¬ 
RQs publicized on the Astronomers TelegrairQ. We did 
not include BLLs in the final analysis because no emis¬ 
sion lines were detected in their optical spectra. Since 
the purpose of the current study is to understand broad 
line variability, we only include objects for which at least 
five epochs of spectroscopy are available for a total of 7 
FSRQs. Tablc[T]lists the SaMOSA sample, Fermi identi¬ 
fier, redshift, number of observations and emission lines 
included in the analysis. 

2 . 2 . The SMARTS Optical/Near-infrared Photometry 
and Optical Spectroscopy 

The SMARTS optical/near-infrared (OIR) photome¬ 
try is obtained nightly with the 1.3m + ANDICAM, 
a dual-channel imager with a dichroic that simultane¬ 
ously feeds an optical CCD and infrared IR imager, with 
spectral coverage from 0.4 - 2.2 /im. Da ta analysis for 
SMAR TS OIR photometry is described in lBonning et alJ 
< 2012 lb so we briefly summarize here. SMARTS differ¬ 
ential photometry is obtained by using optical compari¬ 
son stars calibrated to Landolt standards on photometric 
nights. The reported SMARTS infrared magnitudes are 
calibrated us ing Two Micron All Sky Survey (2MASS) 
magnitudes dSkrutskie et all 12006H of at least one sec¬ 
ondary star in the same field as the blazar. In both op¬ 
tical and infrared bands, the photometric uncertainties 
are dominated by the random errors in the comparison 
star magnitude, and the la uncertainty is reported. In 
the infrared, the dominant uncertainty is the calibration 
to the 2MASS magnitude. OIR finding charts and com¬ 
parison star magnitudes for PKS 0402-362, PKS 0454- 
234, PKS 2142-75 and PKS 2052-474 are provided in 
the Appendix; SMART S finding charts for th e remaining 
sources can be found in lBonning et all ( 2012 1. OIR pho¬ 
tometry and finding charts for all SMARTR-monitored 
sources are publicly available via our websitc@- 

SMARTS optical spectroscopy was obtained with the 
1.5m + Cassegrain spectrograph (RCSpec) at an //7.5 
focus with plate scale 18."1 mm -1 and a LORAL IK 
(1200 x 800) CCD. The primary grating for this study 
has first order resolution of 17.2 A, with spectral cov¬ 
erage of 6600 A and 2 " slit width. In the optimal 
case, spectroscopic data were obtained approximately bi¬ 
weekly, depending on weather conditions and source vis¬ 
ibility in Cerro Tololo, although in many cases signifi¬ 
cantly less data was acquired. No second order correc¬ 
tions were applied to the spectra, as previo us analysis 
yielded a contamination rate of 8 % or less (llsler et al.l 
l2013f) , which is insignificant compared to the errors intro¬ 
duced by the flux calibration. Nevertheless, the system¬ 
atic error introduced by the contaminating continuum 
light is difficult to model, especially in the case when the 
variability amplitude is larger in the blue than the red 
part of the continuum. In any case, this contamination 

1 http://www.astronomerstelegram.org 

2 http://www.astro.yale.edu/smarts/glast/home.php 
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applies only to emission lines with observed wavelength 
> 7000 A, so that Ha in PKS 1510-089 is the only emis¬ 
sion line affected by this contamination. 

The optical spectrosco pic data re d uction here is simi¬ 
lar to that descri bed bvllsler ct al l (120131b We use the 
MPFIT package (lMarkwardtll2009r i to measure emission 
line equivalent width by fitting a Gaussian to the emis¬ 
sion line above the continuum and minimizing the x 2 
statistic. A linear fit to the continuum was applied on 
each side of the line over 100 A. We calculated the noise 
per pixel by combining the uncertainties from bias sub¬ 
traction, sky correction, and aperture extraction. 

The uncertainty in the equivalent width was estimated 
by running 500 Monte Carlo simulations of each fitted 
line, including the measured noise in the count rate of 
each pixel (as described above). For each Monte Carlo 
simulation, the emission line was fitted and the values of 
the equivalent width were calculated. The reported error 
of the equivalent width of each line is the standard error 
on the 500 equivalent width simulations. 

Emission line flux was derived from the equivalent 
width and SMARTS B-, V- or R-band photometry, de¬ 
pending on the observed location of the line in the 
spectrum. Emission lines with rest-frame line cen¬ 
ter, A obs < 5000 A were calibrated with the B-band 
(A e ff = 4400 A), line centers 4999 < A 0 b s < 5999 A 
were calibrated with the V-band (A e ff = 5500 A), 
and X 0 bs > 6000 A were calibrated with the R-band 
( A V7 = 6600 A) . No observed line was more than 500 A 
from the relevant effective wavelength in any case. Ta¬ 
bic [ 2 ] lists the targets, observations, emission line equiv¬ 
alent widths with their associated uncertainties, as well 
as B-, V-, R- and J-band magnitude with associated ler 
uncertainty for the SaMOSA sample. 

2.3. Steward Observatory Optical Spectroscopy and 
Polarimetry 

Since the launch of Fermi, the Steward Observatory 
of the University of Arizona has carried out regular, 
publicly available optical spectrophotometry and lin¬ 
ear spectropolarimetry of a large sample of 7 -ray-bright 
blazars using the Bok 2 .3m and Kuiper 1.54m tele¬ 
sco pes dSmith et apl2009lh The SPOL spectropolarime- 
ter dSchmidt et al.lll992l f is used for this monitoring pro¬ 
gram. Observations are made in first order using a 600 
1 mm -1 grating, yielding a spectral range of 4000 - 7550 
A and resolution of 15-20 A. 

Sky-subtracted spectra and broad-band (5000-7000 A) 
polarization measurements derived from the spectropo¬ 
larimetry from 2008-2013 were obtained from the pub¬ 
licly accessible Steward Observatory website^- The op¬ 
tical spectroscopy was then reduced in a similar fashion 
as the SMARTS spectroscopy. The noise per pixel array 
for the Steward spectroscopy was estimated a posteriori 
using the gain and read noise. This array was then scaled 
by a factor of 10 to approximate the root-mean-square in 
the observed spectrum. The scaling of the uncertainties 
provides a very conservative estimate of the error in the 
lines, and likely overestimates the noise in the spectrum. 
The broad-band polarization measurements were used to 
derive V-band polarized flux densities by multiplying the 

3 http:/www.james.as.arizona.edu/~psmith/Fermi/ 


fractional linear polarization (P) by the V-band flux den¬ 
sity (Fy). 

There is a small systematic offset, of order(O.l) dex, in 
the line fluxes reported by Steward and SMARTS , likely 
due to differences in the comparison stars used for the 
analysis. The (logarithmic) offsets from the SMARTS 
data are applied to the light curves of PKS 0402-362 
(Mg II: -0.1203), PKS 0454-234 (Mg II: -0.2493), and 
PKS 1510-089 (H 7 : +0.0481, H:-0.0501); the other four 
sources included in the SAMOSA sample contain only 
SMARTS data so no cross-calibration is necessary. 

2.4. Fermi 'y-ray Fluxes 

Fermi /LAT data were obtained from the first 
SMARTS photometric observation for each source 
through 2013 July 01 (MJD 56474), via the Fermi Sci¬ 
ence Support Center websitefl Pass 7 reprocessed data 
(event class 3) were analyzed using Fermi Science Tools 
(v9r33p0) with scripts that automate the likelihood anal¬ 
ysis. Galactic diffuse models (gll_iem_v05_revl), isotropic 
diffuse background (iso_p7v6source) and instrument re¬ 
sponse functions (P7REP_CLEAN_V15) were utilized in 
the analysis. Data were constrained to time periods 
where the zenith angle was less than 100° to avoid Earth 
limb contamination, and photons to within a 10 ° region 
centered on the source of interest. 

The Fermi 7 -ray spectra of each object were mod¬ 
eled as a power law or log-parabola, according to spec¬ 
tral type listed in the 2FGL catalog, with the photon 
flux and spectral index as free parameters. Fermi light 
curves were integrated in one-day bins to match the av¬ 
erage SMARTS photometric cadence, and an integral 
Fermi 7 -ray flux above lOOMeV, F 7 is reported. Fermi 
fluxes for which TS > 16 are plotted in subsequent fig¬ 
ures, where TS is the Fermi test statistic and y/TS is 
roughly eq uivalent to the detection significance p er inte¬ 
grated bin dMattox et al.lll996l : IAbdo et al.ll2009 + When 
daily binned fluxes fell below the significance threshold, 
we plot weekly binned flux, also at the TS> 16 level. In 
the case of PKS 0454-234, adaptive binning techniques 
dLott et al.l 120121) were used to determine the Fermi 7 - 
ray flux in the 10 day span around MJD 56280. The same 
TS threshold that was used to plot significant Fermi 7 - 
ray fluxes during the daily and weekly analysis. 

We note that simultaneous measurements of Fermi and 
Steward data means data obtained within the same day. 
The SMARTS optical spectroscopy was flux-calibrated 
to OIR photometry to within the hour. In no case are 
any two datasets matched with temporal separation of 
more than 4 days. 


3. RESULTS 

We evaluate the emission line variability of the 
SaMOSA sample by two independent methods. First, 
we define an empirical emission line flare as a significant 
excursion of the line flux above the mean level, foll ow¬ 
ing a prescription similar to that presented in INalcwaikol 
(|2013f) . but adapted to the present dataset. Specifically, 

4 http://fermi.gsfc.nasa.gov/ssc/data/access 
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Table 1 

SaMOSA Observation Summary 


Source Name 

Fermi Identifier 

Redshift No. of Observations Observed Line(s) 

PKS 0208-512 

2FGL J0210.7-5102 

1.003 

37 

Mg II, C III] 

PKS 0402-362 

2FGL J0403.9-3604 

1.423 

9 

Mg 11, Si III], C IV 

PKS 0454-234 

2FGL J0457.0-2325 

1.003 

43 

Mg 11 

3C 454.3 

2FGL J2253.9+1609 

0.859 

35 

Mg 11, H 7 , H/3, Ha 

PKS 1510-089 

2FGL J1512.8-0906 

0.36 

102 

Mg 11, H 7 , H/3, Hat 

PKS 2052-474 

2FGL J2056.2-4715 

1.489 

8 

Mg 11, C III], C IV 

PKS 2142-75 

2FGL J2147.4-7534 

1.139 

17 

Mg II 


Note. — I An emission line that may have second order contamination. See text for details. 


Table 2 

SaMOSA Optical Photometry and Spectroscopy Log 


UTCt 

MJD 

B 


V 

(T V 

R 

°R 

J 


W(Mg II) 

a W(Mg II) 

W(C III]) 

a w( cm]) 

PKS 0208-512 

20080623 

54640.379 

16.472 

0.004 



15.49 

0.004 



4.953 

0.467 

9.859 

1.672 

20080805 

54683.331 

18.016 

0.015 

17.497 

0.018 

17.137 

0.018 

15.77 

0.023 

- 

- 

13.321 

4.131 

20080823 

54701.3 

18.155 

0.024 

17.654 

0.025 

17.282 

0.022 

- 

- 

18.807 

2.828 

16.451 

3.933 

20080908 

54717.236 

18.372 

0.045 

17.654 

0.058 

17.433 

0.048 

- 

- 

23.452 

3.344 

27.091 

6.037 

20080926 

54735.151 

16.738 

0.007 

16.283 

0.008 

15.882 

0.007 

14.494 

0.011 

6.5 

0.927 

8.273 

1.567 


Note. — tpjTC is in YYYYMMDD format. The equivalent widths of the emission lines are reported as W(species), e.g. W(Mg II), 
in units of Angstroms. SMARTS photometry are given in magnitudes. Table [2] is published in its entirety in The Astrophysical 
Journal, a portion is shown here for guidance regarding its form and content. 



Figure 1 . The daily binned Fermi (E>100MeV) 7 -ray photon 
flux (purple circles), or, when undetected in daily binning, weekly 
binned photon flux (light purple squares), as well as emission 
line fluxes of Mg II (magenta circles) and C III] (cyan stars) for 
PKS 0208-512. 

we define an empirical emission line flare if 1) at least 
three consecutive points, with any two sequential points 
separated by < 60 days, and 2) at least one point be¬ 
tween the first point to deviate from the mean and the 
last point to deviate from the mean must be > 3cr above 
the mean line flux. The empirical line flares are recorded 
in Table |4] 

The second method for evaluating emission line vari¬ 
ability is the y 2 variability test. We compute the y 2 
statistic for an assumed constant line flux (the null hy¬ 
pothesis) and use the number of degrees of freedom 
(DOF) to calculate the probability that the given y 2 de¬ 
viates from the null hypothesis. Unless otherwise noted 



MJD 

Figure 2. Light curves for PKS 0402-362. Symbols as in Fig¬ 
ure [I] plus C IV (green triangles), Si III] (orange stars). Steward 
Observatory data are also presented in the Mg il light curve (open 
pink circles), offset by -0.1203 dex to match the SMARTS mean 
line flux measurement. As the total variability (deviation from the 
respective mean) is evaluated, the normalization does not impact 
the results. The large Fermi 7 -ray flare at approximately MJD 
55830 is not well sampled in the emission line light curve. 

for an individual source, the y 2 variability test was con¬ 
sistent with the null hypothesis at the p> 0.05 level. Re¬ 
sults can be seen in Table [5] This method is sensitive 
to the calculated uncertainty in the emission lines. The 
SMARTS spectra have better-constrained uncertainties 
than the Steward data, which were derived a posteriori , 
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Figure 3. Light curves for PKS 0454-234, with symbols as in 
Figure 1. Optical linear polarimetry from Steward Observatory 
(grey squares), reflects the non-thermal (rather than total) flux 
contribution. A significant emission line flare in Mg II, peaked at 
log F Mg II = -14.23 erg s —1 cm -2 on MJD 56285, is indicated 
by the vertical grey dot-dashed line. Emission line fluxes from 
Steward Observatory are represented by open circles and offset by 
-0.2493 dex. 



2 4 6 8 10 

MJD - 56280 


Figure 4. Fermi 7 -ray light curves for PKS 0454-234, centered on 
the significant Mg II emission line flare. Hourly bins (l avender cir¬ 
cles; top), derived using the adaptive binning technique (ILott et all 
120m The fluxes produced via the adaptive binning method meet 
the same TS threshold as the daily binned fluxes (TS > 16). 


Table 3 

Mean Emission Line and 7 -ray Fluxes 
and Luminosities 


Source 

Line 

(L line ) 

(L/ine) 

(L 7 ) 

PKS 0208-512 

Mg II 

-14.17 (0.06) 

43.94 

51.29 (0.20) 


C III] 

-14.18 (0.18) 

43.95 


PKS 0402-362 

C IV 

-13.90 (0.06) 

44.20 

51.74 (0.28) 


S III] 

-14.10 (0.12) 

44.01 


Mg IT 

-14.34 (0.08) 

43.38 


PKS 0454-234 

Mg IT 

-14.64 (0.07) 

43.08 

51.35 (0.24) 

PKS 1510-089 

Mg II 

-13.51 (0.10) 

43.14 

50.62 (0.31) 


H 7 t 

-13.92 (0.19) 

42.72 



H/3+ 

-13.65 (0.20) 

42.97 



Ha 

-13.20 (0.08) 

43.44 


PKS 2052-474 

Mg 11 

-14.09 (0.09) 

43.48 

51.75 (0.19) 


C IV 

-14.55 (0.20) 

44.07 



C III] 

-14.82 (0.44) 

43.67 


PKS 2142-75 

Mg 11 

-14.18 (0.10) 

43.69 

51.58 (0.20) 

3C 454.3 

Mg 11 

-13.96 (0.10) 

43.60 

52.17 (0.39) 


H/3 

-13.93 (0.14) 

43.61 



H 7 

-13.87 (0.14) 

43.67 



Note. — Mean line flux and line luminosity is given in 
units of log erg s —1 cm -2 and log erg s —1 , respectively, and 
includes all available data. 'The systematic offset between 
the Steward and SMARTS data has been applied before 
calculating the mean line fluxes and luminosities. Fermi 
7 -ray mean luminosity includes all Fermi data with TS > 
16 in the entire observing window for each source and is in 
units of log erg s — b The standard deviation is the reported 
uncertainty. Dat a for 3C 454.3 were previously reported in 
llslcr et all 12013D . but are reproduced here for comparison. 


Table 4 

Significant Emission Line Flares 


Source 

Line 

MJDt 

UTC* 

Line Fluxt 

&line 

s* 

PKS 0454-234 

Mg II 

56285 

20121222 

-14.23 

0.15 

3.0 

PKS 1510-089 

Ha 

54934 

20090411 

-13.09 

0.01 

8.5 


Ha 

55292 

20100404 

-13.12 

0.02 

3.8 


H 7 

55292 

20100404 

-13.79 

0.04 

3.3 


H/3 

55292 

20100404 

-13.58 

0.02 

3.7 


H/3 

55622 

20110430 

-13.57 

0.02 

4.0 


H/3 

56050 

20120501 

-13.55 

0.02 

5.4 

3C 454.3** 

Mg II 

55165 

20091130 

-13.73 

0.06 

1.8 


H 7 

55165 

20091130 

-13.58 

0.09 

2.8 


H 7 

55518 

20101118 

-13.43 

0.06 

3.7 


Note. — Emission line fluxes are in units of log erg s” 1 cm -2 . 
'The MJD and associated emission line flux is given for the peak 
value. tfJTC is given in YYYYMMDD format. ’Detection sig¬ 
nificance, S, is given in units of a away from the mean line flux. 
**The significances for 3C 454.3 are reproduced here for compar¬ 
ison. 


so while any variability detected in the SMARTS data 
using this method was not negated by the addition of 
Steward data, the statistical significance was weakened 
in some cases due to the larger uncertainty on a given 
line flux measurement. For this reason, we calculate the 
X 2 statistics based only on the SMARTS data, for which 
the uncertainties are well-characterized. 

Figures [I][ 8 ] show the Fermi 7 -ray and emission line 
light curves for the SaMOSA sample. 

Notes on Individual Sources 
3.1. PKS 0208-512 

PKS 0208-512 was observed with SMARTS from 2008 
August 6 - 2013 February 12 (MJD 54640 - 56339) and 


the emission line behavior is shown in Figure[T] Although 
a slight increase in line flux was detected in C III] be¬ 
tween MJD 54701 - 54892, peaking at MJD 54735, with 
log F CIII i = -13.97 erg s -1 cm -2 (2.5ct), it does not meet 
the empirical emission line flare significance criteria for 
an empirical emission line flare and is not included in 
subsequent analysis. Thus, no significant emission line 
variability was seen in this blazar. 

3.2. PKS 0+02-362 

PKS 0402-362 was observed from 2011 October 2 - 2012 
December 12 (MJD 55836 - 56285). The source under¬ 
went a large, short-duration 7 -ray flare from MJD 55821 - 
55838 that is not well-sampled in emission line flux (see 
Figure H. No significant emission line variability was 
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Figure 5. Light curves for PKS 1510-089. Symbols as in previous Figures, plus H 7 (orange stars), H/3 (green squares) and Ha (blue 
diamonds). Four well-defined emission line flares are observed, as indicated by the grey dot-dashed lines. In the Fermi 7 -ray flaring period 
from MJD 554850 - 55000, Ha underwent a strong line flare that peaked on MJD 54934 at F Ha = -13.97 erg s —1 cm -2 . During the same 
7 -ray flare, TeV photons were detected by H.E.S.S. between MJD 54910 - 55952. H 7 and H (3 show line flares on MJD 55292 at 3.3cr and 
3.7cr, respectively. 
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MJD 

Figure 6. The light curve for PKS 1510-089, marked as in Figure |5j now shown with the four significant emission line flaring periods 
isolated for each emission line. Dates are given in units of MJD - 55000 in all but the first panel, which are given in MJD - 54000. In each 
case, there is a Fermi 7-ray flare (or increased 7-ray emission) associated with each emission line flare. However, not every emission line 
had detectable emission line variability. 
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MJD 


Figure 7 . Light curves for PKS 2052-474. Symbols as in previous 
Figures. 



log F y [ph s' 1 cm -2 ] 



MJD 


Figure 8. Light curves for PKS 2142-75. Symbols as in previous 
Figures. 


Table 5 

X 2 Variability Statistics 


Source 

Line 

DOFt 

X 2 

Prob. 

cr 

PKS 0208-512 

Mg II 

33 

27.8 

0.72 

0.35 


C III] 

34 

59.6 

0.004 

2.86 

PKS 0402-362 

Mg II 

4 

1.39 

0.84 

0.20 


Si III 

4 

3.13 

0.54 

0.62 


C IV 

4 

5.03 

0.28 

1.07 

PKS 0454-234 

Mg II 

4 

12.9 

0.01 

2.51 

PKS 1510-089 

Mg II 

17 

18.4 

0.36 

0.91 


H 7 

15 

32.3 

0.01 

2.76 


H /3 

17 

87.9 

< 10“ 3 

6.75 


H a 

14 

199.2 

< 10“ 3 

12.3 

PKS 2052-474 

Mg II 

9 

17.0 

0.07 

1.78 


C III] 

11 

7.78 

0.73 

0.34 


C IV 

9 

3.94 

0.91 

0.11 

PKS 2142-75 

Mg II 

11 

14.5 

0.21 

1.26 

3C 454.3 

Mg II 

24 

61.9 

< 10" 3 

4.14 


H /3 

25 

86.9 

< 10- 3 

5.64 


H7 

26 

43.2 

0.02 

2.36 


Note. — t DOF are the degrees of freedom, x 2 
is the total statistic. Prob. is the probability that 
the flux deviates from the null result and cr is the 
significance from the mean, assuming a normal dis¬ 
tribution. 


Figure 9 . Differential optical polarized flux vs. Fermi 7-ray flux 
for PKS 0454-234. The correlation coefficient for the linear regres¬ 
sion is 0.37. Polarization is not a “clean” measure of the strength 
of the non-thermal emission, given its dependence on both the or¬ 
der of the magnetic field strength and the source brightness. For 
this source, the non-thermal synchrotron flux present in the optical 
spectrum is essentially unrelated to the intensity of the 7-ray flux. 



log F t [ph s 'em 2 ] 


Figure 10 . Optical polarized flux vs. Fermi 7-ray flux for 
PKS 1510-089. The correlation coefficient for the linear regression 
is 0.60. The non-thermal synchrotron flux in the optical spectrum 
is not very strongly correlated with the 7-ray flux but a trend is 
clear. 
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Figure 11 . Optical polarized flux vs. Fermi 7 -ray flux for 
3C 454.3. The correlation coefficient for the linear regression is 
0.59. Although the correlation is weak, we see that the polariza¬ 
tion does increase with increasing 7 -ray flux. 

detected in this blazar. 

3.3. PKS 0454-234 

PKS 0454-234 was observed between 2011 October 16 
- 2013 July 27 (MJD 55850 - 56500); the light curve can 
be seen in Figured We also plot the optical linear polar¬ 
ization. For this object, calibrated V-band photometry 
from Steward Observatory was not available, so we used 
differential photometry (AFy) to determine the relative 
polarized flux. Mg II underwent a 3.0cr line flare with 
peak line flux log F Mgii = -14.23 erg s " 1 cm " 2 on MJD 
56285. No accompanying Fermi 7 -ray flare was detected 
in the dai l y binned flux es, so the adaptive binning tech¬ 
nique dLott et al.l I2012T) was utilized during the 10-day 
period around MJD 56280 to determine if sub-day vari¬ 
ability was present. We plot the results of the adaptively 
binned Fermi 7 -ray fluxes along with the daily binned 
7 -ray fluxes and the Mg II emission line flare for com¬ 
parison in Figured The apparent deviation of 7 -ray flux 
derived from adaptive binning is of the same order as the 
textitrms in this region and is thus not significant. 

3.4. 3C 454.3 

We first presented t h e em ission line light curve for 
3C 454.3 in ITsler et all (fjoTjh : one can also be seen in 
iLeon-Tavares et alJ (120 1 il l for this epoch. We extend the 
previous analyses by applying the criterion for emission 
line flares used here; both previously reported emission 
line flares in H 7 and Mg II meet both variability criterion 
and are thus statistically significant. 

3.5. PKS 1510-089 

Optical spectra were obtained for PKS 1510-089 be¬ 
tween 2008 May 17 - 2012 May 3 (MJD 54603 - 56050); 


Table 6 

Fermi 2FGL Variability Index 


Source Name Variabiity Index 

3C 454.3 

14189 

PKS 1510-089 

6405 

PKS 0454-234 

1501 

PKS 0402-362 

1417 

PKS 2142-75 

1162 

PKS 2052-474 

791 

PKS 0208-512 

733 


Note. — The variability in¬ 
dex is obtained fro m the Fermi 
2 yr Source Catalog INolan et al.| 

120m 

the light curve can be seen in Figure [3] We detected 
line flares in Her peaking on MJD 54934 at log Fho= - 
13.09 erg s " 1 cm " 2 with 8.5<r significance. The lead¬ 
ing line flux for this Ha flare was also above the mean 
at MJD 54913 with 3.3cr significance (log Fho= -13.14 
erg s ” 1 cm" 2 ), suggesting that the emission line flare 
may extend past the range observed on the increasing 
side of the 7 -ray flare. During the Fermi 7 -ray flare 
(MJD 54850 - 55000) with peak flux F 7 = -5.08 ph s " 1 
cm " 2 on MJD 54917, H.E.S.S. also detected very high 
energy photons from MJD 55910 - 55952, with highest 
emission, lo g F> n istyv ~ -10.4 ph s " 1 c m " 2 on MJD 
54915 ([H.E.S.S. Collaboration ctaLll201 31. 

H 7 and H/3 show significant line flares on MJD 55292 
at the 3.3cr (log Fh 7 = -13.79 erg s " 1 cm" 2 ) and 3.7cr 
(log Fh/ 3 = -13.58 erg s " 1 cm" 2 ) level, respectively. A 
significant emission line flare was observed in H/3 on 
MJD 55622 with log Fh/?= -13.57 erg s " 1 cm " 2 (4.0cr). 
H/3 has an emission line flare during the same 7 -ray 
flare, peaking on MJD 56050 at log Fh<*= -13.54 erg 
s " 1 cm" 2 (5.4cr). In Figure [ 6 ] we show the regions where 
emission line flares were detected in at least one emission 
line, as described above. 

Statistically significant emission line variability in the 
Ha, H 7 and H/3 emission lines is measured and the x 2 
variability test also indicates variability: Xhc* = 199.4 (14 
degrees of freedom, p<C 10" 3 ), Xh 7 = 32.3 (15 degrees of 
freedom, p=0.01) and Xh/ 3 = 37-9 (17 degrees of freedom, 

p<10" 3 ). 

3.6. PKS 2052-474 

PKS 2052-474 was observed from 2011 May 2 - 2012 
July 11 (MJD 55683 - 56119) and is the only source in 
the sample that did not show any significant Fermi 7 -ray 
flux above log F 7 = -6 ph s " 1 cm " 2 over the epoch of 
observation presented here (see Figured. No significant 
emission line variability was observed in this blazar. 

3.7. PKS 2142-75 

PKS 2142-75 was observed from 2010 May 5 - 2012 
September 15 (MJD 55321 - 56185), seen in Figure [3] 
No significant emission line flares were observed during 
the epoch of observation. 

4. ACCRETION DISK - JET INTERACTION 

The statistically significant emission line variability 
seen in PKS 0454-234, 3C 454.3 and PKS 1510-089, but 
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not in less 7 -ray active sources, suggests that in the most 
active 7 -ray sources the broad line region may be par¬ 
tially photoionized by the jet. In this case, a correlation, 
potentially with lags, between the 7 -ray flux and emis¬ 
sion line flux is expected. We characterize the 7 -ray jet 
activity by using the variabili ty index , as define d in the 
second Fermi catalog (2FGL; iNolan et aLll2012f) . as the 
sum of 2 x log(likelihood) comparison between the flux 
fitted in 24 time segments and a flat light curve over the 
full 2 yr catalog interval. Values greater than 41.64 indi¬ 
cate that there is less than 1 % chance of being a steady 
source. We find that the three blazars in which we identi¬ 
fied emission line variability in this sample also have the 
highest variability index, although all the sources in this 
sample are consistent with 7 -ray variability with high 
certainty. We list the variability indices in Table [ 6 ] 

A discrete correlation analysis would be useful to test 
such a correlation, however, with the small number of 
data points in a given emission line light curve it is not 
valid for this data set. We attempted to quantify the 
significance of the correlation between Fermi 7 -ray flare 
and emission line flare above a random occurrence by 
carrying out a Monte Carlo simulation. We ran 1000 it¬ 
erations of randomized emission line fluxes with respect 
to the date of spectroscopic observation and then applied 
our two tests of emission line variability to each iteration. 
Because the distribution of the randomized emission line 
fluxes is completely dominated by clustering due to the 
seasonal observation schedule and relatively small num¬ 
ber of observations, we were not able to derive a mean¬ 
ingful significance. 

If this explanation of a possible correlation (potentially 
with temporal offset) between the 7 -ray flux and emission 
line flux is correct, then other empirical indicators should 
confirm non-thermal emission in the optical regime dur¬ 
ing emission line flares. To test the presence of such emis¬ 
sion, we consider the optical polarization of PKS 0454- 
234, 3C 454.3 and PKS 1510-089. Polarization data are 
not available for the other sources in the sample, as they 
are too far south of Steward Observatory. 

While optical photometry measures the total emission 
from both the accretion disk and the jet, the optical 
polarization measures the contribut ion of synchrotron 
flux in the optical -ultraviolet regime dRaiteri et al.ll2012l ; 
ISmith et alJ ll994lh Thus, by calculating the polarized 
flux, we can measure the level of non-thermal emission 
at the same time as the 7 -ray and emission line fluxes 
in order to distinguish between the thermal contribution 
and additional emission from the jet. However, the po¬ 
larized flux depends on both the magnetic field ordering 
of the synchrotron region, and how bright it is. These 
two parameters are not often closely related, and as a 
result, the polarized flux is not a “clean” measure of the 
strength of the non-thermal continuum compared to the 
7 -ray flux (which comes only from the jet). Yet, we can 
associate this emission with non-thermal emission given 
the featureless spectrum and rapid variability of the po¬ 
larized flux, which is too rapid to originate from scattered 
accretion disk thermal emission. 

We plot the polarized flux light curve for PKS 0454- 
234 and PKS 1510-089 i n Figures [3] and [5j respectively, 
and refer the reader to llsler et aKl fl2013n for a similar 
plot for 3C 454.3. Visual inspection confirms the gen¬ 
eral agreement of an increase in polarized flux during 


7 -ray flares, although temporal offsets can be seen be¬ 
tween bands. We plot the polarization with respect to 
7 -ray flux of these 3 blazars in Figures l9l [TUI and [TT1 We 
find a correlation coefficient of R = 0.37, R = 0.59 and R 
= 0.60 for PKS 0454-234, 3C 454.3 and PKS 1510-089, 
respectively. Thus, we suggest that there is an increase 
in polarized light during periods of increased 7 -ray ac¬ 
tivity in 3C 454.3 and PKS 1510-089, while we do not 
infer such a relationship in PKS 0454-234 based on these 
data. 

Taken together, the polarization, 7 -ray and line flux 
diagrams provide tentative evidence that, during 7 -ray 
flaring events, an additional population of non-thermal 
ionizing photons produce enough photoionizing flux to 
increase the emission line flux. We infer this from the 
temporal proximity (but not necessarily simultaneity) of 
the emission line fluxes to 7 -ray flares in the light curve 
and the coincidence of optical polarization during peri¬ 
ods of high 7 -ray flaring and emission line fluxes. We 
note that all of the proxies for correlation here are likely 
diluted by offsets in the timescales of increasing flux in 
the different bands. Still, these empirical results suggest 
that there could be a significant source of photoionizing 
flux being produced by a non-thermal jet, and that this 
emission may be contributing to the increase in emission 
line flux that we observe. 

5. DISCUSSION 

5.1. Emission Line Variability 

We find that 4 of the 7 blazars observed in this sample 
show no evidence of statistically significant emission line 
variability. The lack of emission line variability in these 
blazars is consistent with the standard model that the 
accretion disk is the predominant source of photoionizing 
flux. 

However, PKS 0454-234, 3C 454.3 and PKS 1510-089 
all show statistically significant emission line variabil¬ 
ity, with a mean peak emission line flare significance 
of 4cr.We are only aware of a few blazars with pub¬ 
lished simultaneous Fermi 7 -ray and optical emission 
line data. Simultaneous emission line variability stud¬ 
ies in Fermi- monitored blazars have generated mixed 
results. Among a set of similar 7 -ray and optically 
bright, variable quasars, little emission lin e varia bility 
was detected in PKS 1222+216 or 4C 38.51 (ISmith et alJ 
1 20l ib iFarina et all 12012b iRaiteri et all 1201 2l l . It was re¬ 
ported that FSRQ PKS 1222+216 did not have signif¬ 
icant emission lin e variability during recent Fermi 7 - 
ray flaring events ([Smith et ai.ll201 it IFarina et al.ll2012f) , 
however line variability at the 2 . 6 cr level is evident in 
both datasets. While emission line variability at this 
level would not have met the criteria set forth in this 
work, the lines do show some variation based solely on 
continuum variations. Without reanalyzing the data, we 
note that when the H li ne l uminosity from Steward Ob¬ 
servatory reported in IFarina et all ( 2012 1 is compared to 
the line luminosity obtained from TNG in the same work, 
there is a 17a difference in line luminosity. The latter was 
obtained following a MAGIC triggered observation. 

5.2. Correlated Variability 

Emission line variability has been detected in 
PKS 0454-234, 3C 454.3 and PKS 1510-089, which are 
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also the three sources with the highest Fermi variability 
index. This is suggested by the empirical emission line 
flares and y 2 variability test, correlation information in¬ 
cluding possible lags is harder to assess given the limits 
of the data. Thus, we attempt to relate the degree of 
non-thermal jet contribution to the photoionizing flux in 
the broad line region to optical polarization. It has long 
been known that during periods of increased 7 -ray activ¬ 
ity, the thermal contribution to the optical-ultraviolet 
continuum is swa mped by non-thermal emission (e.g., 
ISmith et al .1 1 19941) . We confirm this result with the op¬ 
tical polarization for 3C 454.3 and PKS 1510-089, where 
brighter 7 -ray (and emission line) fluxes correspond to 
highly polarized states, presumably due to jet emission. 
However, we do not find evidence of this behavior in 
PKS 0454-234. The synchrotron peak has been shown 
to be well -correlated from infrared to ultraviolet for FS- 
RQs (e.g. iBonning et aTll2009l) . such that optical varia¬ 
tions indicate similar variability patterns in the ultravi¬ 
olet, where the photoionizing flux peaks. 

The tentative picture that emerges from this work is 
that the most active 7 -ray flaring sources have evidence 
of emission line variability that is not seen in sources with 
less active jets. The correlations between 7 -ray flux and 
emission line flux are likely diluted by the presence of 
lags and/or leads in the peak of either curve. Thus, we 
expect that correlated line variability could be detected 
with higher cadence, multi-epoch optical emission line 
studies of 7 -ray active blazars. 


5.3. Location of the 7 - emitting Region 


Next, we turn our attention to the location of 
the 7 -emitting region, which has been the subject 
of much research. The models fall into two major 
categories: 1 ) near-field 7 -emitting models suggest 

that the bulk of the jet dissipation occurs on sub-pc 
scales, either very deep in th e Broad Line Region 


the broad line region ( 

Bottcheh 20071 Kataoka et al. 

20081 Tavecchio et al.l 

20101 : Ghisellini et al.l 

2010 ; 

Poutancn & Stern 

2010; Tavecchio et all 

2010 ; 

Stern & Poutancn 

20lit Abdo et al.1 201 lib but in 


any case at or be l ow ca nonical distances of 0.1 pc 
(e.g. iPetersonl 119931 l2006fl . and 2) far-held 7 -emitting 
scenarios, which suggest jet dissipation on much larger 


spacial scales (tens of pc) from the c entral source 
dMarscher et al.1 1201 1:: lAgudo et al.l l20lil l.Iorstad et al. i 
HoH). These studies attempt to constrain the location 
of the 7 -emitting region based on SED modeling, 
correlation studies, and/or ultra-short 7 -ray variability. 


More recently, a few coordinated optical spectroscopic 
variability studies, like the one presented here, have 
been undertaken in an effort to ide ntify emissio n line 


variability of v-rav bright blazars 

Smith et al.l 

2009; 

Benitez et al.l 20101: Raiteri et al. 2012 

: ISmith et al. 

2011 : 

Isler et al.ll2013: Leon-Tavares et al.l I2013D. In the 

cases 


where line variability is found, attempts are sought to di¬ 
rectly (and simultaneously) relate the broad line variabil¬ 
ity to jet variability via 7 -ray, mm or other non-thermal 
emission. 

The results of both the indirect and direct studies have 
been mixed, even for the same source and same hare. For 
example, in 3C 454.3, following the 2009 December and 
2010 November haring periods, iTavecchio et al.l (|2010D 


and lAbdo et~ahl (j2011l) suggest the 7 -emitting region of 
3C 454.3 could be located at the outer edges of the broad 
l ine region (r Pm ~ 0.14 pc), using 77 -opacity arguments. 
Ilsler et ah •: 20 1 :il: also suggested a near-held dissipation 
mechanism after observing statistically signiheant emis¬ 
sion line variab ility in both the 2009 and 2010 hares. By 
contrast, ILeon-Tavares et all (1201311 argue in favor of a 
potentially far-held dissipation mechanism in the 2010 
hare, given statistically signiheant emission line variabil¬ 
ity in close temporal proximity to a mm-core ejection. 
They argue that their lack of detectable emission line 
variability during the 2009 hare, in combination with the 
absence of an additional mm-core ejection, suggests that 
emission line variability may be caused by the radio core 
ejections. The two spectroscopic studies come to differ¬ 
ent conclusions likely due to different observation win- 
do ws around the 2009 hare. The observations presented 
bv ILeon-Tavares et al.l (20131 ) did not extend across the 
entire 7 -ray hare, therefore, a simultaneous comparison 
of the peak 7 -ray to e mission line f luxes was not possible. 
The data collected bv Ilsler et al.l (1201311 extended across 
the entire 7 -ray hare,albeit with fewer total observations. 
Thus, the lac k of detected emissio n l ine variability in the 
2009 hare by iLeon-Tavares et all (120131 ) is likely due to 
lack of temporal coverage and not to the absence of line 
variability itself. 

We also co nsider the mm-core ejections seen in 
PKS 1510-089 dMarscher et al.1(201 Of ) with respect to the 
emission line variability reported here. We report emis¬ 
sion line variability in Ha which peaks on MJD 54934. A 
mm-core ejection is not seen in this haring period until 
MJD 54959, suggesting that it is likely not the cause of 
the emission line variability (on the pc scales on which 
the core ejection is observed). To our knowledge, no 
subsequent studies on core ejections in PKS 1510-089 
have been published. Therefore, we are unable to com¬ 
pare our reports of emission line variability. While these 
two examples of non-coincident core ejections at the time 
of emission line variability do not preclude such occur¬ 
rences, we suggest that a core ejection is not a necessary 
condition for emission line variability. However, the near 
temporal proximity of 7 -ray hares to nearly every in¬ 
stance of emission line variability does suggest that jet 
hares (and the attendant increase in photoionizing hux) 
are likely required to produce emission line variability on 
timescales of a few weeks to months. Thus, we consider 
our results to favor a near-held jet dissipation region. 

5.4. Broad Line Region Structure 

While the emission line hux variability is evidence 
of a BLR response, it does not provide a conclusive 
test of its dynamical structure. Whether the jet emit¬ 
ting region is within the canonical BLR or part of an 
outhowing wind (either from the disk or entrained in 
the jet) on larger scales, is still unknown. The broad 
emission line prohles alone are not sufficient to dis- 
tinguish between bro ad line region dynamical models 
dCanriotti et al.l 11980 ). although higher line profile mo¬ 
ments like asymmetry ca n distinguish steady-s tate ver¬ 
sus dynamical theories (jCapriotti et all flDSll) . Thus, 
whether the BLR is a distribution of gas clouds (or fila¬ 
ments) driven primarily by Keplerian velocities, or a disk 
(or jet) wind cannot be constrained by the current line 
variability studies; higher resolution (and cadence) coor- 
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dinated spectroscopic observations are necessary to truly 
constrain the dynamical nature of the BLR in blazars. 

However, if some part of the BLR were located very 
far from the central source as a result of entrainment in 
the jet, we would expect to see a non-variable line core 
with highly variable (blueward) line wings due to the 
high velocities of these outflowing and entrained clouds. 
Conversely, we do not expect the jet emitting region to 
be located deep within the broad line region, as it would 
require a more isotropic (likely unbeamed) jet contribu¬ 
tion that would not be significantly brighter than the ac¬ 
cretion disk in the source frame and hence not produce 
detectable emission line variability. In addition, observa¬ 
tions of very high energy (TeV) emission coincident with 
emission line variability, like that seen in PKS 1510-089, 
are hard to reconcile with the increased probability for 
77 -absorption deep in the rich BLR photon field. 

If, instead, the jet emission were interacting with the 
broad li ne region, a s sug gested in the ‘mirror model’ 
(iGhisellini fc Madaul 1199611 . one could still observe jet- 
augmented photoionization in the broad line emission. 
According to this model, which estimates the BLR as 
a thin shell, the energy density increases dramatically 
at the location of the BLR, where the gas sees strongly 
beamed flux from the jet (see their Figure 2). Two fac¬ 
tors are at play here: 1 ) the jet emission is beamed 
and illuminates only a small fraction of the assumed 
spherical BLR, within a solid angle 7 t/T 2 , though with a 
strongly enhanced flux, and 2 ) the disk radiation is as- 
sume d to be roughly isotr opic. Together with equation 
26 of iGhisellini fc Madaul (11996IL this implies the mea¬ 
sured BLR luminosity, L blr ~ L disk+^p^, where both 
L disk and L j e t are observed luminosities. For anT~ 10, 
the jet contribution to the total (optical-UV) photoion- 
izing luminosity is a factor of ~ 100 less that that of the 
disk. Thus, for the observed variation of a factor ~ 2 in 
BLR flux, the jet photoionizing flux should increase by 
a factor of > 100 . 

For 3 C 454.3 and PKS 1510-089, this constraint is eas¬ 
ily accommodated, as both have shown significant in¬ 
creases (AB > 3 mag) in optical-UV flux during flaring 
events and L jet > L disk, espe cially in high fla ring s tates. 
However, as was shown in dMarscher et all l2010T) . the 
emitting region in a given blazar is quite complicated 
and has been observed to shift from sub-pc to several 
pc scales in a matter of months (|Finke fc Permed 120101 
IGhisellini et all [201311 and thus we do not expect every 
flare to take place within the BLR. 

In the case where the 7 -emitting region is located out¬ 
side the BLR, we would not expect flares in the jet to 
cause line variability. However, lines can vary because 
of variable disk emission; this should happen on longer 
time scales (since the disk varies slowly compared to the 
jet) and need not, in general, be associated with an in¬ 
crease in 7 -ray emission. The latter kind of event may 
have taken place in PKS 0454-234, where we see insuf¬ 
ficient optical-UV jet photoionizing flux combined with 
the lack of a strong 7 -ray flare near the occurrence of the 
emission line flare suggest that this instance of line vari¬ 


ability may be caused by a different physical mechanism. 
In this case, the accretion disk variability could be caused 
by hotspots which can change the sha pe and behavior o f 
the optical sp ectrum as described bv (IR.uan et al .1120 Mil . 
Furthermore, IGhisellini et all (|2010h has already shown 
that L disk ~ L j e t in PKS 0454-234, such that the jet 
does not produce enough photoionizing radiation in the 
optical-UV regime to produce the necessary photoioniza¬ 
tion of the BLR. We also note that the optical polariza¬ 
tion is not well correlated with 7 -ray flux in PKS 0454- 
234 as in 3C 454.3 and PKS 1510-089, suggesting a dif¬ 
ferent physical mechanism causes the line variability in 
PKS 0454-234. 

6 . CONCLUSIONS 

Over 5 years of observations, we find little statisti¬ 
cally significant emission line variability in the 7 sources 
presented here. However, in 3 7 -ray flaring blazars, 
PKS 0454-234, 3C 454.3 and PKS 15i0-089, significant 
emission line variability was detected. We test the opti¬ 
cal polarization as a proxy for non-thermal jet contribu¬ 
tion during periods of emission line variability in these 
sources. We find that in some cases the optical polar¬ 
ization increases with the Fermi 7 -ray flux and emission 
line flares, but the correlation is poor, probably because 
the details of the jet structure affect the polarization sig¬ 
nal. From the 7 -ray flaring, we infer the presence of 
non-thermal photoionizing photons in the system that 
could interact with the broad line region and cause the 
observed emission line increases. 

While we cannot conclusively determine whether there 
are lags between the emission line increases and the 7 - 
ray flares, due to poor temporal sampling in the present 
dataset, we find that the most 7 -ray active blazars have 
statistically significant line variability that is not seen in 
less 7 -ray active sources. 

Higher cadence optical spectroscopy is needed to inves¬ 
tigate better the correlation between the emission lines 
and the 7 -ray fluxes. This requires near daily spec¬ 
troscopy in both active and quiescent states to build up 
enough data to truly constrain the degree of correlated 
variability between the jet and emission line flux. 
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Figure 12 . The SMARTS optical (left) and infrared (right) finding charts for PKS 0402-362, PKS 0454-234, PKS 2052-474, and PKS 2142- 
75. The field of view is 5.12' X 5.12' for the optical finding charts and 1.45' X 1.45' for the infrared. Comparison star magnitudes and 
lcr uncertainties are given in Table [7] optical comparison stars are labeled with numbers and infrared comparison stars are labeled with 
letters. In cases where optical and infrared comparison stars are the same, numbers are used to identify the star in both images. 




Figure 13 . Optical and infrared finding charts for PKS 0454-234, labels as in Figure [121 

APPENDIX 

SMARTS OPTICAL AND INFRARED FINDING CHARTS 

The SMARTS optical and infrared finding charts for PKS 0402-362, PKS 0454-234, PKS 2052-474, and PKS 2142-75 
are presented here with the BVRJK magnitudes used to calibrate the com parison stars. SMART S OIR finding charts 
for PKS 1510-089 and PKS 0208-512 have been previously published in IBonning et al.l (120121) . All SMARTS OIR 
finding charts can be found on our website. 
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Figure 14 . Optical and infrared finding charts for PKS 2052-474, labels as in Figure WI\ 




Figure 15 . Optical and infrared finding charts for PKS 2142-474, labels as in Figure fl2l 
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Table 7 

Optical and Infrared Comparison Stars 


Target 

Star 

B (<r B ) 

V (cr v ) 

R (o-fl) 

J (crj) 

K (cr K ) 

PKS 

0402-362 

i 

15.38 

(0.03) 

14.75 

(0.03) 

14.34 

(0.02) 

13.42 

(0.03) 

13.02 

(0.03) 



2 

18.47 

(0.06) 

17.14 

(0.04) 

16.18 

(0.02) 

- 


- 




3 

19.17 

(0.09) 

17.81 

(0.03) 

16.81 

(0.03) 

- 


- 




A 

- 


- 


- 


16.82 

(0.14) 

15.82 

(...) 

PICS 

0454-234 

1 

18.04 

(0.04) 

17.02 

(0.03) 

16.32 

(0.03) 

- 


- 




2 

17.53 

(0.03) 

16.55 

(0.02) 

15.86 

(0.02) 

- 


- 




3 

18.12 

(0.04) 

17.05 

(0.03) 

16.32 

(0.02) 

- 


- 




4 

19.22 

(0.06) 

17.83 

(0.03) 

16.82 

(0.03) 

- 


- 




5 

16.65 

(0.02) 

16.00 

(0.02) 

15.62 

(0.02) 

- 


- 




A 

- 


- 


- 


15.60 

(0.06) 

14.86 

(0.12) 



B 

- 


- 


- 


13.45 

(0.02) 

12.65 

(0.03) 

PKS 

2052-474 

1 

18.29 

(0.05) 

17.87 

Tonry 

17.52 

(0.05) 

- 


- 




2 

17.51 

(0.02) 

16.94 

(0.01) 

16.59 

(0.03) 

- 


- 




3 

15.28 

(0.02) 

14.67 

(0.01) 

14.29 

(0.02) 

- 


- 




4 

16.43 

(0.02) 

15.82 

(0.01) 

15.46 

(0.03) 

- 


- 




5 

16.99 

(0.02) 

16.23 

(0.02) 

15.80 

(0.03) 

- 


- 




6 

16.25 

(0.02) 

15.56 

(0.01) 

15.15 

(0.01) 

- 


- 




A 

- 


- 


- 


11.44 

(0.02) 

11.24 

(0.02) 



B 

- 


- 


- 


15.78 

(0.11) 

14.43 

(0.10) 



C 

- 


- 


- 


16.43 

(0.11) 

16.67 

(...) 

PKS 

2142-75 

1 

17.25 

(0.04) 

16.12 

(0.02) 

15.41 

(0.02) 

13.90 

(0.03) 

12.96 

(0.04) 



2 

19.11 

(0.08) 

18.28 

(0.03) 

17.75 

(0.04) 

16.51 

(0.14) 

15.84 

(...) 



3 

17.23 

(0.04) 

16.51 

(0.02) 

16.04 

(0.03) 

- 


- 




4 

19.23 

(0.06) 

18.39 

(0.04) 

17.90 

(0.04) 

- 


- 




5 

17.98 

(0.04) 

16.93 

(0.03) 

16.18 

(0.03) 

- 


" 



Note. — The reported uncertainties are la. Optical comparison star magnitudes are 
listed by number, and infrared comparison star magnitudes are listed by letter. When 
the same comparison star is used in the optical and infrared, the requisite data is labeled 
by number for both bands. Infrared uncertainties in 2MASS that were not available via 
the catalog are denoted by These data are calibrated to the 23 June 2014 optical 
and infrared photometry. 
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